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Abstract 

We have performed extensive molecular dynamics simulations to study noise power spectra of 
density and potential energy fluctuations of a Lennard-Jones model of a fluid in the supercritical 
region. Emanating from the liquid-vapor critical point, there is a locus of isobaric specific heat 
maxima, called the Widom line, which is often regarded as an extension of the liquid-vapor coexis- 
tence line. Our simulation results show that the noise power spectrum of the density fluctuations 
on the Widom line of the liquid- vapor transition exhibits three distinct 1//'^ behaviors with ex- 
ponents 7= 0, 1.2 and 2, depending on the frequency /. We find that the intermediate frequency 
region with an exponent 7 ~ 1 appears as the temperature approaches the Widom temperature 
from above or below. On the other hand, we do not find three distinct regions of 1//'^ in the 
power spectrum of the potential energy fluctuations on the Widom line. Furthermore, we find that 
the power spectra of both the density and potential energy fluctuations at low frequency have a 
maximum on the Widom line, suggesting that the noise power can provide an alternative signature 
of the Widom line. 

PACS numbers: 61.20.Ja, 64.70.F-, 05.40.Ca, 64.60.Bd 
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I. INTRODUCTION 



In a typical pressure-temperature (P — T) phase diagram for a fluid, there is a first 
order phase transition line between the liquid and the vapor phases that terminates at a 
critical point [1]. Beyond this critical point lies the supercritical regime where one can go 
continuously between the liquid and vapor phases. In this region thermodynamic response 
functions such as the isobaric specific heat Cp, the isothermal compressibility kt and the 
thermal expansion coefficient a do not monotonically increase or decrease with P and T. 
Rather they exhibit a line of maxima {2 4| emanating from the critical point that is known 
as the Widom line Q which can be thought of as an extension of the liquid-vapor phase 
boundary. 

More generally, a Widom line is regarded as an extension of a first order phase boundary 
that terminates at a critical point. It is marked by maxima in the thermodynamic response 
functions. The most studied Widom line is in supercooled water where a first order liquid- 
liquid transition has been proposed between a high-density liquid and a low-density liquid 



phase 



-|8|. This phase boundary has been hypothesized to end in a critical point from 



which the Widom line emanates jSKZl]. This critical point cannot be directly investigated 
experimentally because it is usurped by the spontaneous crystallization of water. Exper- 
imental searches jol (foj for this critical point have involved putting water in a confined 



ll|-|l4|. However, so far there has been no 



geometry to avoid spontaneous crystallization 
direct experimental evidence of this critical point. Therefore, studies focused on the Widom 
line have been based on the belief that a Widom line is associated with the existence of a 
critical point {2, O]- These studies have found a rich and complex behavior. For example, 
along the Widom line in supercooled water, experiments have found a fragile-to-strong dy- 



namic transition 



151], a sharp change in the proton chemical shift which is a measure 



of the local order in confined water 



16|, and a density minimum of water jol. Simulations 



of supercooled water have found a breakdown of the Stokes-Einstein relation [ITI and a 

n n 

fragile-to-strong dynamic crossover [5|, Il8| along the Widom line. Simulations of biological 
macromolecules with hydration shells find that the glass transition of the macromolecules 
coincides with water molecules in the hydration shell being on the Widom line [l9 |. 

Given the interesting behavior associated with the liquid-liquid phase transition of su- 
percooled water, we decided to focus on the liquid- vapor Widom line which has been largely 



neglected. Experiments have found that the hquid-vapor Widom hne significantly affects 
dynamic quantities as well as static quantities. X-ray diffraction measurements of the struc- 
ture factor S{Q) of fluid argon have found a liquid-like phase with high atomic correlation as 
the Widom line is approached with decreasing pressure to something intermediate between 
a liquid and a gas with low atomic correlation 20|. (The amount of atomic correlation is 
indicated by the first sharp diffraction peak in S{Q) versus Q.) Inelastic X-ray scattering 
experiments on the sound velocity in fluid oxygen find liquid-like behavior above the Widom 
line 2l|, and, for the case of fluid argon, a sound velocity that increases with wavevector 



at pressures in the liquid-like phase above the Widom line but not below it 



22|. Recent 



theoretical simulations of a Lennard- Jones (LJ) fluid and a van der Waals fluid have found 
maxima in Cp, k^, a, and density fluctuations along the Widom line [s], [4, |23 |. 

In this paper we use noise spectra to study the liquid-vapor Widom line. Noise spectra 
have been used to probe first and second order phase transitions where it has been found 
that the low frequency noise is a maximum at the phase transition 2J|. In addition, in the 
vicinity of second order transitions the noise goes as l/f^ (/: frequency) where the exponent 



7 can be related to the critical exponents |24|-|27l|. Motivated by this, we decided to probe 



28| . Away from 



the supercritical region and the Widom line with the noise power spectra of fluctuations 
in the density and the total potential energy generated from 3D molecular dynamics (MD) 
simulations of an LJ fluid. Along the Widom line we find that the noise power spectrum Sp 
of the density p goes as l/f^ where the exponent 7 = (white noise) at low frequencies, 
7 = 1.2 at intermediate frequencies, and 7 = 2 at high frequencies. The behavior at 
intermediate frequencies is indicative of a distribution of relaxation times 
the Widom line, there is no intermediate regime, only white noise at low frequencies and 
7 = 2 at high frequencies. At low frequencies where there is white noise in both the density 
and the potential energy, the magnitude of the low frequency noise is a maximum along 
the Widom line. This is similar to the maximum in the low frequency noise in the energy 
and magnetization that was found at the second order phase transition temperature of the 
2D Ising ferromagnetic Ising model and at the first order transition of the 2D 5-state Potts 



model 



24| . A maximum in the low frequency noise is consistent with the maxima of the 



thermodynamic response functions at the phase transitions and along the Widom line. Here 
we present MD simulation studies of the relation between the Widom line and the power 
spectra of the density and potential energy fluctuations of the supercritical LJ fluid. 



This paper is organized as follows. In Sec. [TTl we describe our MD simulation of an LJ 
fluid in detail. In Sec. IIIII we present the results of our calculations of thermodynamic 
quantities. We show the pressure-temperature and density-temperature phase diagrams, 
and we calculate the thermodynamic response functions, such as the thermal expansion 
coefficient and the isobaric specific heat, to find the location of the Widom line in the phase 
diagram. We also calculate the standard deviations of the density and potential energy 
distributions, and find that they have their maximum value along the Widom line. In 
Sec. IIVI we present the noise power spectra of the density and potential energy fiuctuations 
of the supercritical L J fiuid. We find 1// noise in the density fiuctuations in the intermediate 
frequency range around the Widom line in addition to white noise at low frequency and 1//^ 
noise at high frequency. We also show that the power spectra of the density and potential 
energy fiuctuations at low frequency have a maximum along the Widom line. In Sec. |V] we 
summarize our results. 



II. METHODS: MOLECULAR DYNAMICS SIMULATIONS OF A LENNARD- 
JONES FLUID 



The LJ interaction has been commonly used in simulations to model fiuids and in par- 
ticular, to study the liquid- vapor phase transition 29|-|33| . The LJ interparticle potential 



has a combination of a short-ranged repulsive core and a long-ranged attractive tail: 



4e 



12 



(1) 



Here the LJ parameters e and a represent the energy and length scales, respectively, and r 
is the distance between two particles. Monte Carlo simulation studies 



30 



31( of LJ fiuids 



have shown that the calculation of the liquid-vapor critical point is sensitive to how the LJ 
interaction is truncated at large distances. For the full LJ interaction without truncation, 
the estimated liquid-vapor critical point is at T* ~ 1.316 and p* ~ 0.304 in reduced units 
(see below) |3o|. For the truncated and shifted LJ interaction, the estimated liquid- vapor 



30 



3l|. 



critical point is located at T* ~ 1.085 and p* ~ 0.317 

We truncated the LJ interparticle potential at a cutoff radius Vc- In order to have the 
interparticle potential and force be continuous at r = Tc, we truncated and then shifted the 



LJ interparticle potential 34| : 



f/(r) 



f/Lj(r)-t/Lj(rc) 



dr 

for r > r, 



(r — Tc) for r < Tc 

(2) 



To reduce the deviation of U{r) from the original LJ interparticle potential, we used a larger 
cutoff radius Vc = 4. Oct than the usual cutoff radius Vc = 2.5a [34 1. 

Throughout this work we use reduced MD dimensionless units with the energy, length 
and mass scales set by e, a, and mass m, respectively. Thus, the reduced length r* = r/a, 



energy E* = E/e, time t* = t/ y^ma'^/e, temperature T* = ksT/e, volume V* = V/a^, 
number density p* = pa^, pressure P* = Pa^/e, specific heat c*p = cp/k^ and thermal 
expansion coefficient a* = ae/k-Q where k-Q is Boltzmann's constant sj]. 

To study the Widom line and power spectra of supercritical LJ fluids, we performed 
extensive MD simulations of an LJ fluid in ensembles with NPT kept fixed where N = 1728 
particles, P is the pressure and T is the temperature. The simulations were done at 33 
temperatures T* from LOO up to L64 with temperature increments of AT* = 0.02 and for 
40 pressures P* from 0.01 up to 0.4 with pressure increments of AP* = 0.01. In the initial 
configuration the particles are arranged in a cubic lattice, then the temperature is raised to 
a very high temperature to melt the lattice, and then the system is brought to the desired 
temperature, and 2 x 10^ MD steps (or equivalently, t* = 10^) used to equilibrate the system, 
followed by 4 X 10^ MD steps (or equivalently, t* = 2 x 10"^) to collect equilibrium data. The 
number of steps need to achieve equilibrium was determined by seeing how many steps were 
needed for the mean potential energy to become time independent. Particle motions were 
updated at each time step At* = 0.005. We implemented the Berendsen thermostat and 
barostat to keep the temperature and pressure constant 35^. We used periodic boundary 
conditions in the x, y and z directions. Since the system is self- averaging, we only needed 
to do one run for any given set of conditions. 



III. RESULTS: THERMODYNAMICS 

The P* — T* phase diagram is shown in Figure [U The first order phase boundary 
between the liquid and vapor phases terminates at the liquid-vapor critical point {T* ~ 
1.305 and P* ^ 0.16). The Widom line continues the liquid-vapor phase boundary beyond 
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the critical point. In the subcritical region where there is a first order phase transition, 
the density p* along the isobaric path abruptly changes when T* crosses the transition 
temperature (Fig. However, p* in the supercritical region continuously chang es as a 
function of T*. These different behaviors of p* as a function of T* are clear when we 
calculate the absolute density difference between two adjacent temperatures |Ap*(T/)| = 
\p*{T*) — p*{T*_^)\ (Figs. |3]and|l]). Note that this differs from the order parameter of the 
liquid- vapor phase transition Ap = puquid — Pvapor- When there is a discontinuous drop in 
p*(T*) at the transition temperature in the subcritical region, |Ap*| shows a spike when T* 
crosses the same transition temperature. In the supercritical region, as shown in Fig. HJ 
|Ap*| shows a maximum at a given temperature, whereas p* continuously changes without 
any abrupt jump. 

The maximum in |Ap*| is associated with a maximum in dp/dT that, as we will see 
later, is associated with the Widom line. Similar behavior has been found in supercooled 
water. As water is cooled, the density reaches a well known maximum at 4 °C. As the 
temperature continues to decrease into the supercooled regime, the density decreases until 



it reaches a minimum as revealed by recent inelastic neutron scattering experiments {q]. 
With further cooling the density increases in the same way that simple liquids do upon 
cooling [36:]. Between the density maximum and minimum, there is a maximum in dp/dT 



which coincides with the Widom line of the liquid-liquid transition {9]. Therefore, as we 
will see later, we expect that the temperature of the maximum in |Ap*| corresponds to the 
Widom temperature T*. 

n 

To find the Widom line, we calculate the thermal expansion coefficient a defined as |1[ 



1 fdV\ 



91np 

dT 



p 

(3) 



Generally, a is associated with fluctuations in the entropy S and volume V 36| . 



< {5S5V) >= VkBTa. (4) 

In Figure El we show a* as a function of T* along the isobaric path in the supercritical region. 
At low temperature, a* initially increases gradually upon heating, and then starts to increase 
rapidly upon further heating. As T* increases further, a* finally reaches maximum, and then 
rapidly decreases. As P* increases above P* (for P* > P*), the magnitude of the peak in a* 
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decreases, and the temperature of the peak in a* moves toward higher T* as P* increases. 
a* has the same temperature behavior as |Ap*| at constant P* . In particular, the maximum 
in a* occurs at the same temperature as that of |Ap*|. 

Next, we calculate another thermodynamic response function, the isobaric specific heat 
cp which is defined as [l| 

1 (dH\ 



p 



1 (d{E + PV)\ 



N 



p 

where H is the enthalpy and E is the internal energy. Generally, cp is associated with 
fluctuations in the entropy S 36| . 



< [dSy >= NkBCp. (6) 

In Figure El we show c*p as a function of T* along the isobaric path in the supercritical region. 
Cp also shows the same temperature behavior that we found for |Ap*| and a* at constant 
P* as shown in Figs. H] and [51 As P* increases from P* (for P* > P*), the magnitude of 
the peak in c*p decreases, and the temperature of the peak in c*p moves toward higher T* as 

;he Cp maximum in Fig. [6] is the Widom temperature T* 



P* increases. The temperature o: 
at a given pressure 



emperc 

m 
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'unctions, the 
IQ. Based 



38| . Note that of all the various response 
isobaric specific heat is usually used to define the location of the Widom line 
on our calculation of Cp, we plot the estimated locations of Cp maxima (the Widom line) in 
Fig. [H denoted by open circles. 

Next we investigate the standard deviation of an observable X defined as 



[{x-xy), (7) 



where X is the average of X. We calculate the standard deviations S^* and Hp* of the 
potential energy per particle u* and density p*, respectively. In Figure [TJ we present Hu* 
and Sp. as a function of T*. The fluctuations in u* and p* rapidly increase as T* approaches 
from below at fixed P*, and they reach their maximum around T* ~ T^. After crossing 
T^, S^j* and Sp* rapidly decrease as shown in Fig. [71 The location in the phase diagram of the 
maxima in the variances agrees well with the location of the maxima of the thermodynamic 
response functions along the Widom line. This is consistent with Eqs. (jlD and ([6|) which 
relate fluctuations to thermodynamic response functions. 



IV. NOISE POWER SPECTRAL ANALYSIS 



Let us set up our notation and define what we mean by noise. Let m{t) be a quantity that 
fluctuates in time. Let 6m{t) be the deviation from its average value of some quantity m at 
time t. If the processes producing the fluctuations are stationary in time, i.e., translationally 
invariant in time, then the autocorrelation function of the fluctuations {5m{t2)5m{ti)) will 
be a function ip{t2 —ti) of the time difference. In this case the Wiener-Khint chine theorem 
can be used to relate the noise spectral density Sm{oj) to the Fourier transform ipiu) of 



the autocorrelation function 



39( 1 : S'm(w) = 2%l)m{!^) where u is the angular frequency. In 



practice Sm{oj) typically is calculated by multiplying the time series 5m{t) by a windowing 
or envelope function so that the time series goes smoothly to zero, Fourier transforming the 
result, taking the modulus squared, and multiplying by two to obtain the noise power 40|. 
(We flnd that our results are not sensitive to the choice of windowing function, so we do not 
use a windowing function; this is equivalent to a rectangular window.) 

1// noise, where / ( = — ) is frequency, corresponds to Smijjj) ~ 1/u. It dominates 
at low frequencies and has been observed in a wide variety of systems, such as granular 
systems, molecular liquids, ionic liquids, a lattice gas model, and resistors |2J, |28|, I4l|-l50|. 



For example, long time fluctuations of the potential energy in water and silica exhibiting 
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48|. The 1// power 



1// spectra have been reported in computer simulation studies |44j . 
spectrum of the potential energy fluctuations for water is related to the non-exponential 
relaxation of slow hydrogen bond dynamics Q,|43] . A study of solvation dynamics in an ionic 
liquid at room temperature has also shown 1// spectral behavior 50|. Fluctuations in the 



number of defects in a disordered two-dimensional liquid also exhibit a 1// power spectrum 



46, 



49| 



at low temperatures, suggesting that the dynamics of the system is heterogeneous 
In addition to the relation between 1// spectral behavior and dynamics, the power spectra 
can be used as a probe of phase transitions. It has been shown that at a phase transition in 
classical spin systems (such as the Ising model and Potts model), the low frequency noise of 
the energy and magnetization fluctuations has a maximum at the transition temperature 24 1. 



A simple way to obtain 1/f noise was given by Dutta and Horn 28|]. We can use the 
relaxation time approximation to write the equation of motion for 5p: 



d6p 
It 



6p 

T 



(8) 
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where r ^ is the relaxation rate. The autocorrelation function V'p(^) is given by 



The Fourier transform is a Lorentzian if there is just one value of r. 

1 + W^T^ 



(9) 



(10) 



where A is an overall scale factor. If there is a broad distribution g{T) of relaxation times, 
then the Fourier transform is a sum of Lorentzians: 



Spiu) = A 



1 + u;V2 



g{r)dT 



(11) 



where T2 and ri are the upper and lower limits of the distribution g{T) and, hence, of the 
integral over r. If we assume that (^(r) ~ r^^ [28], then we obtain 

1 A f I 



Sp{uj) oc — 

CO 



arctan ( — arctan ( . (12) 

In the low frequency limit (/ ^ Tg"^ ^ ''"f^)) the power spectrum shows 1/ behavior, and 
in the high frequency limit (/ ^ ^ ^ 'T-z'^)^ shows 1//^ behavior. In the intermediate 
frequency region (r2"^ <^ f <^ '^1^)1 the power spectrum exhibits 1// behavior 28|. 

Here we use the noise power spectra of p* and probe of the Widom line. The 

power spectrum of an observable X{t) is defined as 

2 



Sxir 



(13) 



There is no subtraction of the mean X, so there will be a delta function at / = 0. Sxif) is 
otherwise the same as the noise spectrum of the fluctuations 6X. 



A. Noise Spectra of Density Fluctuations 

In Figures [8] to [T3|, we present the power spectra Sp* of the density fluctuations as a 
function of frequency /* for different temperatures along the isobaric path of P*= 0.20. 
When the temperature T* (= 1.02) is far below the Widom temperature ~ 1.36 such 
that T* ^ T^, Sp* shows two different frequency dependences: almost flat (~ 1//*°) at low 
frequencies and rapidly decreasing (~ 1//*^) at high frequencies (see Fig. |8]). At T* = 1.26, 
the thermodynamic response functions start to change rapidly, and Sp* (/*) clearly exhibits 
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three different frequency dependences: flat (~ 1//*°) at low frequencies, relatively slow 
decreasing (~ 1//*^'^) with increasing frequency at intermediate frequencies, and rapidly 
decreasing (~ 1//*^) at high frequencies (see Fig. On the Widom line (T* = 1.36 ~ T^), 
Sp* , as shown in Fig. [TUl shows the same three different frequency dependences as in Fig. [HI 
flat (~ 1//*°) at low frequencies, relatively slow decreasing (~ l/f*^'"^) at intermediate 
frequencies and rapidly decreasing (~ 1//*^) at high frequencies. For different P*, Sp*{f) ~ 
1//*'^ at the corresponding Widom temperature T^{P*) with 7= 0, 1.2 and 2, respectively 
(see Fig. [TTj) . Interestingly, in the intermediate frequency region, 7 ~ 1.2 all along the 
Widom line, independent of P*. 



The behavior of Sp* at temperatures in the vicinity of is similar to that seen along the 
Widom line. Below but close to where the thermodynamic response functions still change 
rapidly (as in the case of T* = 1.44 shown in Fig. [T^ . Sp* shows three different frequency 
dependences: white (~ 1//*°) at low frequencies, relatively slow decreasing (~ 1//*^ '') at 
intermediate frequencies and rapidly decreasing (~ 1//*^) at high frequencies. When T* is 
far away from T^, there are two frequency dependences: flat (~ 1//*°) at low frequencies 
and rapidly decreasing (~ 1//*^) at high frequencies (see Fig. [T3] where T* = 1.58 ^ ^w))- 



When T* is far away from T^, the density fluctuations relax exponentially with a single 
relaxation time. As explained in the simple model of 1// noise in Eqs. ffTOj) and f|T2|) . 1/f 
spectral behavior can be connected with a broad distribution of relaxation times. Therefore, 
our results imply that as T* approaches T^, the density fluctuations are becoming more 
complicated and are relaxing with a broad distribution of relaxation times. The appearance 
of the 1/f spectral behavior near the Widom line in a narrow frequency range might be 
from the fact that the Widom line of the liquid-vapor transition is located at very high 
temperatures, at which heterogeneous dynamics generally does not occur 46|, |49| , so that 
the distribution of relaxation times might have a narrow range of the values as implied by 
Eq. f ll2p . Note that Sp* along the Widom line does not provide a clear distinction between 



liquid-like and vapor-li 
of sound in argon did 



ffi behaviors as the recent x-ray scattering experiments on the velocity 
21. 
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B. Noise Spectra of Potential Energy Fluctuations 



We now examine the fluctuations in the potential energy u* per particle. In contrast 
to Sp*, the noise spectra 5*^. of potential energy fluctuations does not clearly show three 
different frequency dependences on the Widom line, as shown in Fig. [T31 Su*{f) is white at 
low frequencies and then slowly decreases with increasing frequency. Previous simulations of 
argon at T = 95 K found that the noise spectrum of the total potential energy fluctuations is 
white over a wide range of frequencies, followed by a rapid decrease at high frequencies 44 1. 
Interestingly, as opposed to the total potential energy, the potential energy fluctuations 
of an individual argon atom exhibit 1// noise in an intermediate frequency range (1 < 
/ < 10 cm"^) [44]. The potential energy of an argon cluster also exhibits 1// noise in an 
intermediate frequency range (/ ^ 1 cm~^), indicating that locally the system can 



nave a 



distribution of relaxation times, e.g., at the core and surface of the argon cluster 4^. The 



origin of the difference between the total potential energy fluctuations of simple liquids and 
the potential energy fluctuations individual atoms and clusters remains elusive 4^. An 
investigation of this difference would be interesting, but is beyond the scope of the present 
work. We should also mention that a previous computational study found that the long-time 
total potential energy fluctuations of water does exhibit 1// noise, unlike simple liquids like 



argon 



44l |. Water has a random hydrogen-bond network that has non-exponential relaxation 



processes, whereas simple liquids do not have such structures 



C. Block Average of Variances in Density and Potential Energy 



We are able to estimate the correlation time of an observable X{t) by calculating the 



block average 



24 



31 



34 



52 



-|54||. To calculate the block average, we take a time series of 



X(t), divide the time series into blocks or bins of equal size, calculate the thermodynamic 
quantity such as the variance of X for each bin, and then average over all the bins to 
obtain the block average. If the quantity is not a linear function of X, then the block average 
can depend on the bin size. So one can use the same time series over again to calculate the 
block average with a different bin size, and see if the block average changes with bin size. 
Typically, the block average will increase with bin size before saturating at a constant value 
corresponding to the equilibrium value of the quantity [2J, |5J]. The saturation occurs when 



11 



the bin size exceeds the correlation time 



34 



52 



53|. 



In Figure [T51 we calculate the block average of the variance in the density (S^,) and 
in the potential energy per particle (S^, ) as a function of the block size (At*) for different 
temperatures. By comparing the block averages for three different temperatures, we find that 
the correlation times of both p* and u* are largest at the Widom temperature (T^ ~ 1.36), 
where |Ap*| and the specific heat have their maxima. The block averages for both quantities 
saturate at a constant value at approximately At* ~ 10^. This block size indicates the largest 
correlation time and corresponds to the crossover frequency between the 1//*° and 1//* 
spectral behaviors that occurs at approximately /* ~ 10~'^. Recent Monte Carlo studies of 
the classical spin systems also found that the crossover between white noise at low frequencies 
and 1/f noise at higher frequencies occurred at a frequency corresponding to the inverse of 
the largest correlation time 2J]. It is interesting to note that at the Widom temperature, 
the correlation time of the density fluctuations is larger than that of the potential energy 
fluctuations. 



D. Maximum of Low Frequency Noise as a Signature of the Widom Line 

A maximum in the low frequency white noise can be a signature of a phase transition 
or a crossover. For example, Monte Carlo simulations have found that the low frequency 
white noise in the energy and magnetization is a maximum at the phase transition of the 2D 



ferromagnetic Ising model and the 2D 5-state Potts model 2J]. In addition, the maximum 



in the noise power of the defect density at low frequency occurs at the onset of defect 



proliferation in a two-dimensional liquid modeled by a Yukawa potential |46l . |49 |. 

We have found that the low frequency white noise in the density and the potential energy 
is a maximum at the Widom line. In Figure [161 "we show values of the power spectra at 
low frequency. We take the average value of the white noise part of the power spectrum Sx 
where X is p* or u* at each T*. Each Sx{T*) value in Fig. [11] is normalized by Sx{Tq = 1) 
to obtain the ratio Sx(T*)/ Sx(T^ = 1). For both the density (X = p*) and potential energy 
(X = u*) fluctuations, Sx(T*)/ Sx{Tq = 1) has a maximum in the vicinity of the Widom 
temperature T^. Therefore, our results show that a maximum in the low frequency power 
can provide an additional signature of the Widom line. 
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V. DISCUSSION AND SUMMARY 



We have performed molecular dynamics simulations of a Lennard- Jones fluid in the su- 
percritical region of the phase diagram to probe the Widom line of the liquid-vapor phase 
transition and its effect on the power spectra of the density and potential energy fluctuations. 



Extending from the liquid-vapor critical point, the Widom line, the locus of the maxima of 
the thermodynamic response functions, is a continuation of the liquid-vapor phase transition 
line. The thermodynamic response functions in the supercritical region, such as the thermal 
expansion coefficient and the isobaric specific heat, have maxima along the Widom line. 
Similar results have been found by studies of the Widom line of the hypothesized liquid- 



liquid phase transition in supercooled water 



m 



19 



37 



38 



55 



56|. 



We studied the power spectra of the density and potential energy fiuctuations in the 
supercritical region. Far away from the Widom line, the noise in the density fiuctuations is 
white at low frequencies and goes as 1//^ at higher frequencies. 1//^ behavior is consistent 
with an exponential relaxation process characterized by a single relaxation time. In the 
vicinity of the Widom line we found that the density noise spectrum could be divided into 
3 frequency regimes: 1//° at low frequency, 1//^'^ at intermediate frequencies, and 1//^ at 
high frequencies. The intermediate region 1/f noise implies that there is a distribution of 
relaxation times associated with the maxima in the response functions along the Widom line. 
Eq. f|T2|) implies that the narrowness of the frequency range where there is 1// noise is due 
to a narrow distribution of relaxation times. The narrow width of the intermediate region 
might be due to the fact that density correlations are short lived at high temperatures. In 
contrast to the density fiuctuations, the power spectrum of the potential energy fiuctuations 
along the Widom line does not exhibit three distinct frequency regimes. Finally, we found 
that the low frequency white noise of the density and potential energy fiuctuations have 
their maxima along the Widom line. This suggests that noise power spectra, which have 
been used to probe phase transitions {2^, can also be used to locate the Widom line. 
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FIG. 1. (Color online) Pressure-temperature (P* — T*) phase diagram of the liquid, vapor and 
supercritical fluid phases. The liquid and vapor phases are denoted by black squares and red 
triangles, respectively. The supercritical fluid phase is denoted by blue crosses. The liquid-vapor 
critical point estimated from our simulation results and Ref. [sO] is denoted by a green diamond. 
The approximate location of the liquid- vapor critical point in the P* — T* diagram is T* ~ 1.305 
and P* ~ 0.16. Circles represent the Widom line emanating from the liquid-vapor critical point. 
The Widom line is the locus of isobaric specific heat maxima Cp estimated from Fig. [6l The Widom 
line is a continuous extension of the liquid-vapor phase boundary. All points here represent phase 
points we simulated. 
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FIG. 2. (Color online) Isobaric trajectories in the density-temperature {p* —T*) phase diagram from 
pressure P* = 0.01 up to P* = 0.29 in increments of AP* = 0.02. A green diamond represents 
the estimated liquid-vapor critical point. Below the critical point, when the system crosses the 
transition temperature, p* shows a discontinuous jump at the transition temperature. Above 
the critical point, p* shows the continuous change over the entire temperature range that we 
investigated. 
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FIG. 3. (Color online) The absolute value of the density difference |A/9*| as a function of temper- 
ature T* for different pressures in the subcritical and supercritical regions. The absolute value of 
the density difference is defined by \Ap*(T*)\ = \p*{T*) — p*{T*_j^)\, which is different from the 
definition of the typical order parameter of the liquid- vapor phase transition Ap* = pJiquid — /Ovapor- 
The approximate location of the liquid- vapor critical point is T* ~ 1.305 and P* ~ 0.16. Here we 
used the temperature step AT.* (= T* — T*_^) = 0.02. 
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FIG. 4. (Color online) Absolute value of the density difference |Ap*| versus temperature T* for 
different pressures in the supercritical region. The magnitude of the peaks decreases as the system 
moves away from the liquid- vapor critical point. At higher P*, the location of the peak moves 
toward higher T* . 
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FIG. 5. (Color online) Thermal expansion coefficient a* versus temperature T* for different pres- 
sures in the supercritical region. As P* increases above the critical pressure P* , the magnitude of 
the peak in a* decreases. At higher P* , the location of the peak moves toward higher T* . 



21 




FIG. 6. (Color online) Isobaric specific heat Cp versus temperature T* for different pressures in 
the supercritical region. As P* increases above P*, the magnitude of the peak in Cp decreases. At 
higher P*, the location of the peak moves toward higher T*. The behavior of c*p is the same as 
that of I A/9* I and a* (see Figs. [Hand ED. 
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FIG. 7. (Color online) Standard deviations (a) S^* of the potential energy u* per particle and (b) 
Tip* of density p* as a function of temperature T*. Both S^* and Sp* show maxima around the 
Widom line. 
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FIG. 8. (Color online) Noise power spectrum Sp* of the density p* as a function of frequency /* 
at T* = 1.02 (< T*) and P* = 0.20 in a log-log plot. The two sohd lines are fits to 1//*'^ with 
the exponents 7= and 2, respectively. 



24 



10"^ 10"^ 10"^ 10"^ 10° 

* 

f 

FIG. 9. (Color online) The noise power spectrum Sp* of the density p* versus frequency /* at T* 
= 1.26 (< r^) and P* = 0.20 on a log-log plot. The point (T* = 1.26, P* = 0.20) is located in the 
liquid region of the phase diagram where the response functions (a and cp) change rapidly below 
Tw. That is, a and cp at T* = 1.26 rapidly increase as T* increases, as shown in Figs. [5]and[6l 
The three solid lines are fits to 1/ f*^ with the exponents 7= 0, 1.1, and 2. 
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FIG. 10. (Color online) Power spectrum Sp* of the density p* as a function of frequency /* at T* 
= 1.36 (~ r^) and P* = 0.20 on a log-log scale. The three solid lines indicate fits to 1//*'*' with 
the exponents 7= 0, 1.2, and 2, respectively. 
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FIG. 11. (Color online) Shown on a log- log scale are power spectra Sp* of the density p* as a 
function of frequency /* at (a) T*= 1.34 and P*= 0.18, (b) T*= 1.38 and P*= 0.22, and (c) T*= 
1.40 and P*= 0.24. These temperatures are the Widom temperature T^{P*) at a given pressure 
P* . Sp* along the Widom line exhibits three distinct frequency regimes where l/f^ has exponents 
7= 0, 1.2 and 2. 
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FIG. 12. (Color online) The power spectrum Sp* of the density p* as a function of frequency /* 
at T* = 1.44 (> r„) and P* = 0.20 on a log-log scale. T* = 1.44 is located in the region of the 
phase diagram where the response functions (a and cp) change rapidly above T^. That is, a and 
cp at T* = 1.44 rapidly decrease as T* increases, as shown in Figs. [5]and[6l The three solid lines 
are fits to 1//*t with the exponents 7 = 0, 1.4 and 2, respectively. 
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FIG. 13. (Color online) Shown on a log-log plot is the power spectrum Sp* of the density p* as a 
function of frequency /* at T* = 1.58 (» T*) and P* = 0.20. The two sohd hues are fits to l/f*^ 
with the exponents 7= and 2, respectively. 
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FIG. 14. (Color online) The power spectrum Su* of the potential energy per particle u* as a 
function of frequency /* at T* = 1.34 (~ T^) and P* = 0.20 on a log-log scale. The two solid lines 
are guides to the eye and indicate 1//*t with the exponents 7= 1 and 2, respectively. 



30 



CM 



5x10 



4x10 



3x10 



2x10 



1x10 



6x10 



5x10 



4x10 



3x10 



2x10 



1x10 




FIG. 15. (Color online) Block averages of (a) the density variance S^, , and (b) the potential energy- 
variance as a function of the block size At* at P* = 0.20 for different temperatures T* = 1.02 

(< r;), 1.36 (~ r;), and i.ss (> r;). 
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FIG. 16. (Color online) The magnitude of the normalized power spectra SxiT*)/SxiTQ = 1) 
versus temperature T* for the density {X = p*) and potential energy per particle {X = u*) at low 
frequency at P* = 0.2. Both Sp*{T*)/Sp*{T^ = 1) and Su- (T*) / Su* {T^ = 1) exhibit maxima in 
the vicinity of the Widom line. 
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